(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     16356,        437]
NotebookOptionsPosition[     15803,        415]
NotebookOutlinePosition[     16213,        431]
CellTagsIndexPosition[     16170,        428]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{"Clear", "[", 
  RowBox[{
  "M", ",", "\[Eta]", ",", "\[Theta]", ",", "\[Lambda]", ",", "\[Gamma]", ",",
    "f"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"fLSO", "[", "M_", "]"}], ":=", 
    FractionBox["1", 
     RowBox[{
      SuperscriptBox["6", 
       RowBox[{"3", "/", "2"}]], "\[Pi]", " ", "M"}]]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Constants", "\[IndentingNewLine]", "\[Gamma]"}], " ", "=", " ", 
     "0.5772"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Lambda]", " ", "=", " ", 
     RowBox[{
      RowBox[{"-", "1987"}], "/", "3080"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Theta]", " ", "=", " ", 
     RowBox[{
      RowBox[{"-", "11831"}], "/", "9240"}]}], ";"}], "\[IndentingNewLine]", 
   "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "0"], "=", " ", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "1"], "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "2"], "=", " ", 
   RowBox[{
    FractionBox["20", "9"], 
    RowBox[{"(", 
     RowBox[{
      FractionBox["743", "336"], "+", " ", 
      RowBox[{
       FractionBox["11", "4"], "\[Eta]"}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "3"], "=", 
   RowBox[{
    RowBox[{"-", "16"}], " ", "\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "4"], "=", " ", 
   RowBox[{"10", 
    RowBox[{"(", 
     RowBox[{
      FractionBox["3058673", "1016064"], "+", 
      RowBox[{
       FractionBox["5429", "1008"], "\[Eta]"}], " ", "+", " ", 
      RowBox[{
       FractionBox["617", "144"], 
       SuperscriptBox["\[Eta]", "2"]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "5"], "=", 
   RowBox[{"\[Pi]", 
    RowBox[{"(", 
     RowBox[{
      FractionBox["38645", "756"], "+", 
      RowBox[{
       FractionBox["38645", "756"], 
       RowBox[{"Log", "[", 
        FractionBox["f", 
         RowBox[{"fLSO", "[", "M", "]"}]], "]"}]}], "-", 
      RowBox[{
       FractionBox["65", "9"], "\[Eta]", 
       RowBox[{"(", 
        RowBox[{"1", " ", "+", " ", 
         RowBox[{"Log", "[", 
          FractionBox["f", 
           RowBox[{"fLSO", "[", "M", "]"}]], "]"}]}], ")"}]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "6"], "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      FractionBox["11583231236531", "4694215680"], "-", 
      FractionBox[
       RowBox[{"640", " ", 
        SuperscriptBox["\[Pi]", "2"]}], "3"], "-", 
      FractionBox[
       RowBox[{"6848", " ", "\[Gamma]"}], "21"]}], ")"}], " ", "+", " ", 
    RowBox[{"\[Eta]", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", 
        FractionBox["15335597827", "3048192"]}], "+", " ", 
       FractionBox[
        RowBox[{"2255", " ", 
         SuperscriptBox["\[Pi]", "2"]}], "12"], "-", " ", 
       FractionBox[
        RowBox[{"1760", " ", "\[Theta]"}], "3"], "+", " ", 
       FractionBox[
        RowBox[{"12320", " ", "\[Lambda]"}], "9"]}], ")"}]}], " ", "+", " ", 
    RowBox[{
     FractionBox["76055", "1728"], 
     SuperscriptBox["\[Eta]", "2"]}], "-", 
    RowBox[{
     FractionBox["127825", "1296"], 
     SuperscriptBox["\[Eta]", "3"]}], "-", 
    RowBox[{
     FractionBox["6848", "21"], 
     RowBox[{"Log", "[", 
      RowBox[{"4", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"\[Pi]", " ", "M", " ", "f"}], ")"}], 
        RowBox[{"1", "/", "3"}]]}], "]"}]}]}]}], ";"}], "\[IndentingNewLine]",
 
 RowBox[{
  RowBox[{
   RowBox[{
    SubscriptBox["a", "7"], "=", 
    RowBox[{"\[Pi]", 
     RowBox[{"(", 
      RowBox[{
       FractionBox["77096675", "254016"], "+", 
       FractionBox["378515", "1512"], "-", " ", 
       RowBox[{
        FractionBox["74045", "756"], 
        SuperscriptBox["\[Eta]", "2"]}]}], ")"}]}]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"PsiF", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{"2", "\[Pi]", " ", "f", " ", 
     SubscriptBox["t", "c"]}], "-", " ", 
    FractionBox["\[Pi]", "4"], "+", " ", 
    RowBox[{
     FractionBox["3", 
      RowBox[{"128", " ", "\[Eta]", " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"\[Pi]", " ", "M", " ", "f"}], ")"}], 
        RowBox[{"5", "/", "3"}]]}]], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"k", "=", "0"}], "7"], 
      RowBox[{
       SubscriptBox["a", "k"], 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"\[Pi]", " ", "M", " ", "f"}], ")"}], 
        RowBox[{"k", "/", "3"}]]}]}]}]}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"CForm", "[", "PsiF", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "0"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "1"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "2"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "3"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "4"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "5"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "6"}], "]"}], "]"}], "\n", 
 RowBox[{"Unset", "[", 
  RowBox[{"Subscript", "[", 
   RowBox[{"a", ",", "7"}], "]"}], "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.4731366348830347`*^9, 3.4731367121913958`*^9}, {
   3.473136748737252*^9, 3.473136787368086*^9}, {3.473136903513839*^9, 
   3.473137386316791*^9}, {3.473182274042816*^9, 3.473182327104337*^9}, {
   3.473182376338027*^9, 3.473182696548211*^9}, {3.473184527088834*^9, 
   3.473184634849889*^9}, {3.473186295786938*^9, 3.473186421347254*^9}, {
   3.473186463704578*^9, 3.473186501974931*^9}, {3.473186551444905*^9, 
   3.473186617793223*^9}, {3.473186657497642*^9, 3.473186673146118*^9}, {
   3.473186738070846*^9, 3.473186749070288*^9}, {3.473186785224772*^9, 
   3.473186785933391*^9}, {3.473187326387244*^9, 3.473187337330909*^9}, {
   3.473187376606365*^9, 3.47318737756042*^9}, {3.473187427585463*^9, 
   3.473187429722158*^9}, {3.473187532658927*^9, 3.473187628156488*^9}, {
   3.473187977916741*^9, 3.473188088313969*^9}, {3.473188165150548*^9, 
   3.473188171584812*^9}, {3.47318820723533*^9, 3.473188211881933*^9}, {
   3.473188350213722*^9, 3.473188355023288*^9}, {3.473188392959709*^9, 
   3.4731883936982*^9}, {3.473188472815013*^9, 3.473188482898269*^9}, {
   3.473188731882007*^9, 3.473188758179309*^9}, {3.473188827512216*^9, 
   3.473188892078215*^9}, {3.473188923627317*^9, 3.473189039340044*^9}, {
   3.473189070986659*^9, 3.473189198895409*^9}, {3.473190417598798*^9, 
   3.473190494857902*^9}, 3.47319054166984*^9, {3.473190588004275*^9, 
   3.473190592343626*^9}, {3.473190622388777*^9, 3.473190623232133*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox["\[Pi]", "4"]}], "+", 
  RowBox[{
   FractionBox["1", 
    RowBox[{"128", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"f", " ", "M"}], ")"}], 
      RowBox[{"5", "/", "3"}]], " ", 
     SuperscriptBox["\[Pi]", 
      RowBox[{"5", "/", "3"}]], " ", "\[Eta]"}]], 
   RowBox[{"3", " ", 
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"16", " ", "f", " ", "M", " ", 
       SuperscriptBox["\[Pi]", "2"]}], "+", 
      RowBox[{
       FractionBox["5", "756"], " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"f", " ", "M"}], ")"}], 
        RowBox[{"2", "/", "3"}]], " ", 
       SuperscriptBox["\[Pi]", 
        RowBox[{"2", "/", "3"}]], " ", 
       RowBox[{"(", 
        RowBox[{"743", "+", 
         RowBox[{"924", " ", "\[Eta]"}]}], ")"}]}], "+", 
      FractionBox[
       RowBox[{"5", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"f", " ", "M"}], ")"}], 
         RowBox[{"4", "/", "3"}]], " ", 
        SuperscriptBox["\[Pi]", 
         RowBox[{"4", "/", "3"}]], " ", 
        RowBox[{"(", 
         RowBox[{"3058673", "+", 
          RowBox[{"5472432", " ", "\[Eta]"}], "+", 
          RowBox[{"4353552", " ", 
           SuperscriptBox["\[Eta]", "2"]}]}], ")"}]}], "508032"], "-", 
      FractionBox[
       RowBox[{"5", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"f", " ", "M"}], ")"}], 
         RowBox[{"7", "/", "3"}]], " ", 
        SuperscriptBox["\[Pi]", 
         RowBox[{"10", "/", "3"}]], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "28137439"}], "+", 
          RowBox[{"4975824", " ", 
           SuperscriptBox["\[Eta]", "2"]}]}], ")"}]}], "254016"], "+", 
      RowBox[{
       SuperscriptBox["f", "2"], " ", 
       SuperscriptBox["M", "2"], " ", 
       SuperscriptBox["\[Pi]", "2"], " ", 
       RowBox[{"(", 
        RowBox[{
         FractionBox["11583231236531", "4694215680"], "-", 
         FractionBox[
          RowBox[{"6848", " ", "\[Gamma]"}], "21"], "-", 
         FractionBox[
          RowBox[{"15335597827", " ", "\[Eta]"}], "3048192"], "+", 
         FractionBox[
          RowBox[{"76055", " ", 
           SuperscriptBox["\[Eta]", "2"]}], "1728"], "-", 
         FractionBox[
          RowBox[{"127825", " ", 
           SuperscriptBox["\[Eta]", "3"]}], "1296"], "+", 
         RowBox[{
          FractionBox["5", "12"], " ", 
          SuperscriptBox["\[Pi]", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "512"}], "+", 
            RowBox[{"451", " ", "\[Eta]"}]}], ")"}]}], "-", 
         FractionBox[
          RowBox[{"1760", " ", "\[Eta]", " ", "\[Theta]"}], "3"], "+", 
         FractionBox[
          RowBox[{"12320", " ", "\[Eta]", " ", "\[Lambda]"}], "9"], "-", 
         RowBox[{
          FractionBox["6848", "21"], " ", 
          RowBox[{"Log", "[", 
           RowBox[{"4", " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"f", " ", "M"}], ")"}], 
             RowBox[{"1", "/", "3"}]], " ", 
            SuperscriptBox["\[Pi]", 
             RowBox[{"1", "/", "3"}]]}], "]"}]}]}], ")"}]}], "-", 
      RowBox[{
       FractionBox["5", "756"], " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"f", " ", "M"}], ")"}], 
        RowBox[{"5", "/", "3"}]], " ", 
       SuperscriptBox["\[Pi]", 
        RowBox[{"8", "/", "3"}]], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "7729"}], "+", 
         RowBox[{"1092", " ", "\[Eta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Log", "[", 
          RowBox[{"6", " ", 
           SqrtBox["6"], " ", "f", " ", "M", " ", "\[Pi]"}], "]"}]}], 
        ")"}]}]}], ")"}]}]}], "+", 
  RowBox[{"2", " ", "f", " ", "\[Pi]", " ", 
   SubscriptBox["t", "c"]}]}]], "Output",
 CellChangeTimes->{{3.473136642607893*^9, 3.473136712616763*^9}, {
   3.473136749479863*^9, 3.473136788121644*^9}, 3.473136899284795*^9, 
   3.4731369397479486`*^9, {3.4731369873436117`*^9, 3.473137065705452*^9}, {
   3.473137179772154*^9, 3.473137196248891*^9}, 3.473137287700967*^9, 
   3.473137369388214*^9, 3.473184624598536*^9, {3.473186396854215*^9, 
   3.473186422271174*^9}, {3.473186473591981*^9, 3.473186503461184*^9}, {
   3.473186567506354*^9, 3.473186583382244*^9}, 3.473186618278175*^9, 
   3.473186673853734*^9, {3.473186740102957*^9, 3.473186750317783*^9}, 
   3.473186937170466*^9, 3.473187337641407*^9, 3.473187377933883*^9, 
   3.473187430612142*^9, 3.473187538891963*^9, {3.47318758486434*^9, 
   3.473187629407714*^9}, 3.473187978646819*^9, 3.473188090382256*^9, 
   3.473188172263854*^9, 3.473188212682146*^9, 3.473188356488079*^9, 
   3.473188483452412*^9, {3.473188750468955*^9, 3.473188759104503*^9}, {
   3.473188832606613*^9, 3.473188892830621*^9}, {3.473188939070646*^9, 
   3.473189006457305*^9}, 3.473189041324551*^9, 3.473189088754709*^9, 
   3.473189149649285*^9, 3.473189201139501*^9, 3.473190425947508*^9, {
   3.47319046144891*^9, 3.473190496197964*^9}, 3.47319054829539*^9, {
   3.473190657155944*^9, 3.473190659204326*^9}}],

Cell["\<\
-Pi/4. + (3*(1 - 16*f*M*Power(Pi,2) + \
(5*Power(f*M,0.6666666666666666)*Power(Pi,0.6666666666666666)*(743 + 924*\
\[Eta]))/756. + 
        (5*Power(f*M,1.3333333333333333)*Power(Pi,1.3333333333333333)*(\
3058673 + 5472432*\[Eta] + 4353552*Power(\[Eta],2)))/508032. - 
        (5*Power(f*M,2.3333333333333335)*Power(Pi,3.3333333333333335)*(-\
28137439 + 4975824*Power(\[Eta],2)))/254016. + 
        Power(f,2)*Power(M,2)*Power(Pi,2)*(2467.5541189728633 - \
(6848*\[Gamma])/21. - (15335597827*\[Eta])/3.048192e6 + \
(76055*Power(\[Eta],2))/1728. - 
           (127825*Power(\[Eta],3))/1296. + (5*Power(Pi,2)*(-512 + \
451*\[Eta]))/12. - (1760*\[Eta]*\[Theta])/3. + (12320*\[Eta]*\[Lambda])/9. - 
           (6848*Log(4*Power(f*M,0.3333333333333333)*Power(Pi,0.\
3333333333333333)))/21.) - 
        (5*Power(f*M,1.6666666666666667)*Power(Pi,2.6666666666666665)*(-7729 \
+ 1092*\[Eta])*(1 + Log(6*Sqrt(6)*f*M*Pi)))/756.))/
    (128.*Power(f*M,1.6666666666666667)*Power(Pi,1.6666666666666667)*\[Eta]) \
+ 2*f*Pi*Subscript(t,c)\
\>", "Output",
 CellChangeTimes->{{3.473136642607893*^9, 3.473136712616763*^9}, {
   3.473136749479863*^9, 3.473136788121644*^9}, 3.473136899284795*^9, 
   3.4731369397479486`*^9, {3.4731369873436117`*^9, 3.473137065705452*^9}, {
   3.473137179772154*^9, 3.473137196248891*^9}, 3.473137287700967*^9, 
   3.473137369388214*^9, 3.473184624598536*^9, {3.473186396854215*^9, 
   3.473186422271174*^9}, {3.473186473591981*^9, 3.473186503461184*^9}, {
   3.473186567506354*^9, 3.473186583382244*^9}, 3.473186618278175*^9, 
   3.473186673853734*^9, {3.473186740102957*^9, 3.473186750317783*^9}, 
   3.473186937170466*^9, 3.473187337641407*^9, 3.473187377933883*^9, 
   3.473187430612142*^9, 3.473187538891963*^9, {3.47318758486434*^9, 
   3.473187629407714*^9}, 3.473187978646819*^9, 3.473188090382256*^9, 
   3.473188172263854*^9, 3.473188212682146*^9, 3.473188356488079*^9, 
   3.473188483452412*^9, {3.473188750468955*^9, 3.473188759104503*^9}, {
   3.473188832606613*^9, 3.473188892830621*^9}, {3.473188939070646*^9, 
   3.473189006457305*^9}, 3.473189041324551*^9, 3.473189088754709*^9, 
   3.473189149649285*^9, 3.473189201139501*^9, 3.473190425947508*^9, {
   3.47319046144891*^9, 3.473190496197964*^9}, 3.47319054829539*^9, {
   3.473190657155944*^9, 3.473190659385076*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"D", "[", 
   RowBox[{
    RowBox[{"Exp", "[", 
     RowBox[{"I", " ", "x"}], "]"}], ",", "x"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.473191326561844*^9, 3.473191362713919*^9}}],

Cell[BoxData[
 RowBox[{"\[ImaginaryI]", " ", 
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{"\[ImaginaryI]", " ", "x"}]]}]], "Output",
 CellChangeTimes->{{3.473191346841762*^9, 3.473191363210103*^9}}]
}, Open  ]]
},
WindowSize->{1237, 972},
WindowMargins->{{Automatic, 102}, {15, Automatic}},
FrontEndVersion->"6.0 for Linux x86 (32-bit) (June 19, 2007)",
StyleDefinitions->FrontEnd`FileName[{"Creative"}, "NaturalColor.nb", 
  CharacterEncoding -> "iso8859-1"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 7255, 201, 876, "Input"],
Cell[7848, 226, 5147, 130, 193, "Output"],
Cell[12998, 358, 2309, 37, 175, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15344, 400, 238, 6, 65, "Input"],
Cell[15585, 408, 202, 4, 43, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
